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The Collins effect connects transverse quark spin with a measurable azimuthal asymmetry in the 
yield of hadronic fragments around the quark's momentum vector. Using two different reconstruction 
methods we measure statistically significant azimuthal asymmetries for charged pion pairs in e'^e^ 
annihilation at center-of-mass energies of 10.52 GeV and 10.58 GeV, which can be attributed to the 
fragmentation of primordial quarks with transverse spin components. The measurement was per- 
formed using a dataset of 547 fb^^ collected by the Belle detector at KERB improving the statistics of 
the previously published results by nearly a factor of 20. 



PACS numbers: 13.88.+e,13.66.-a,14.65.-q,14.20.-c 

I. INTRODUCTION 

Electron-positron annihilation at high energies leads 
to the creation of quark-antiquark pairs with high mo- 
menta. Assuming that neither quark in a given event ra- 
diates an energetic gluon, the quark and the antiquark 
fragment into two jets of hadrons in opposite hemi- 
spheres. The energies and the momentum vectors of the 
quarks can be inferred from the energies and momenta 
of the observed final state hadrons. The quark-hadron 
fragmentation process is usually parameterized with the 
help of fragmentation functions D'^{z). These describe 
number densities for a quark of flavor q to fragment into 
a hadron h carrying the fraction z of the original quark 
momentum. 

The two final state quarks can be created in transverse 
spin states. Here we use azimuthal correlations between 
pairs of hadrons in opposite jet hemispheres to study 
transverse spin effects in the quark fragmentation pro- 
cess. Specifically, we study the distribution of final state 
pions around the momentum vector of the fragmenting 
quark. The quark momentum direction is measured ap- 
proximately by identifying it with the thrust axis. 

Spin-dependent effects in the fragmentation of quarks 
into hadrons were first discussed by Collins the 
Collins function H^{z, k^) is the amplitude of the mod- 
ulation in the azimuthal distribution of the final state 
hadrons in spin-dependent fragmentation processes. It 
depends on the normalized hadron momentum z and the 
magnitude of the transverse hadron momentum kx with 
respect to the three momentum of the quark. Initial ef- 
forts to extract Collins asymmetries in e+e^ annihilation 
were carried out in DELPHI • The first observation of 
the azimuthal Collins asymmetries in quark fragmenta- 
tion was reported by the Belle collaboration . This pa- 
per is an extension of the earlier Belle measurement with 
nearly a factor of 20 higher statistics allowing a more re- 
fined analysis. 

The Collins function cannot be calculated reliably em- 
ploying QCD based algorithms, because the fragmen- 
tation process is a non-perturbative process. However, 



once determined experimentally the Collins effect can be 
used as "quark spin polarimeter" to determine the trans- 
verse spin states of final state quarks in semi-inclusive 
deep inelastic scattering (SIDIS) and in polarized proton- 
proton scattering. Over the past few years there has been 
increasing interest in spin phenomena and experiments 
have published new data on the transverse spin struc- 
ture of the nucleon 0, d, IE S IE Hi ■ However, the 
semi-inclusive observables in SIDIS are not a direct mea- 
surement of the transverse quark spin distributions (the 
so-called transversity distributions). Instead they mea- 
sure a product of quark transversity distributions and 
the Collins fragmentation functions. Only an indepen- 
dent measurement of the Collins asymmetry in e+e" an- 
nihilation makes it possible to extract transversity quark 
distributions from single transverse spin asjnnmetries in 
SIDIS or polarized proton-proton scattering. 

Recently, the first global transversity analysis of SIDIS 
and Belle data has been carried out by Anselmino and 
collaborators ifiotl . The resulting u and d quark transver- 
sity distributions have large errors that are dominated 
by the experimental uncertainties in the SIDIS and Belle 
data. The new results presented here in combination 
with new data expected in the near future from HER- 
MES and COMPASS will reduce the uncertainties in the 
extraction of quark transversity distributions. 

The study of transverse momentum dependent (TMD) 
fragmentation and distribution functions such as the 
Collins fragmentation function has led to important the- 
oretical developments in QCD and have greatly advanced 
the understanding of factorization and process depen- 
dence of distribution- and fragmentation-functions in 
hard scattering processes. TMD hard quark scattering 
can give rise to the large single spin asymmetries that 
have been observed experimentally but seemed to be 
incompatible with the traditional hard scattering QCD 
framework, which assumes factorization ifTTll . Specifi- 
cally, factorization assumes that the exchange of soft glu- 
ons in the initial or final state of a hard scattering process 
can be neglected. Furthermore it is assumed that parton 
distributions and fragmentation functions are process- 
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independent, i.e. that the CoUins function in SIDIS will 
be the same as in e+e^. It has now become possible to 
correctly include the initial and final state interactions 
in QCD calculations for some hard scattering processes. 
Recent calculations ifl^ ITsIl with the correct treatment 
of initial and final state interactions in hard processes 
show that the Sivers effect lIT^ll arises from this inclu- 
sion. Another consequence of this new treatment of fac- 
torization is the process dependence of the Sivers distri- 
bution function. While the Sivers distributions in SIDIS 
and the Drell-Yan process are expected to have the same 
magnitude they are predicted to have opposite signs. 

Collins fragmentation appears to be an important test 
case for future theoretical approaches that describe low 
energy QCD phenomena. The Collins function describes 
phenomenologically the evolution into hadrons starting 
from a quark with a given momentum and spin orienta- 
tion. The challenge for the future will be to calculate this 
process from basic principles. 

The paper is organized as follows: after a discussion 
of the Collins effect the experimental procedure will be 
described in detail, a detector description is followed by 
a description of the analysis and the study of systematic 
errors. Finally, results will be presented together with 
an attempt to interpret the measured Collins asymme- 
tries for di-hadrons in terms of the Collins fragmentation 
functions. 



A. Transversity 

Recently there has been an increased interest in the 
third leading twist quark distribution function of the nu- 
cleon, the so-called transversity distribution Sq{x) [fisll . 
The corresponding charge is the tensor charge, which 
is obtained by summing the charge-squared weighted 
transverse spin distribution over all flavors, integrated 
over the quarks' momentum fraction x. Due to its sym- 
metry properties this is a more valence-like object than 
the axial charge (which is obtained in a similar way from 
the quark helicity distribution) [510. Since the transver- 
sity distribution function is chiral-odd, gluons cannot 
contribute at leading order. Therefore no gluon transver- 
sity distribution exists and thus the QCD evolution of 
transversity is different from that of the quark helicities. 
Transversity cannot be accessed with inclusive DIS exper- 
iments, since the corresponding forward scattering am- 
plitude contains a helicity flip that is suppressed by the 
nearly vanishing quark mass relative to the scale \/Q^. 
The scale in DIS is defined as = —q^^q^ where is the 
four-momentum transfer from the lepton to the nucleon. 

As a consequence one can only measure transversity if 
a second chiral-odd distribution or fragmentation func- 
tion is involved. One possibility to measure transversity 
uses Drell-Yan processes, via double spin asymmetries 
of two transversely polarized proton and (anti) proton 
beams. This jdelds information on the product of a 
quark and an antiquark transversity distribution. Un- 



fortunately the cross section for this process is quite 
small and the possibility of reasonably high transverse 
polarizations for antiprotons is limited although some 
efforts are being made to conduct such experiments in 
the future lfl6l IitII . Another way to access transversity 
is by combining it with a chiral-odd quark fragmentation 
function. This combination is available in semi-inclusive 
DIS off a transversely polarized nucleon by detecting at 
least one of the produced hadrons. Such measurements 
are currently being performed at HERMES [4] on a 
proton and at COMPASS |@] on a deuteron target. In 
addition, a very similar observable can be accessed 
in proton-proton collisions where one of the beams is 
transversely polarized. In both cases there are different 
candidates for chiral-odd fragmentation functions, which 
are briefly discussed in the following three paragraphs. 



1. Collins fragmentation function 

Most relevant to the interpretation of existing data on 
experiments with transversely polarized targets or beams 
is the Collins function first proposed in Ref. 1 1 1. The 
Collins function relates the transverse spin of the quark 
to the transverse momentum of a final state hadron. 
These correlations result in a non-uniform azimuthal dis- 
tribution of final state hadrons around the initial quark 
direction. 



2. Interference fragmentation function 

The interference fragmentation function ifTsIl mea- 
sures the interference of two hadrons in the final state 
between two different angular momentum states. This 
interference is again visible as an azimuthal asymmetry 
of the plane defined by the two hadrons around the orig- 
inal quark momentum. According to model predictions 
ifTgl [20I1 the interference fragmentation function may be 
different from zero whenever the invariant mass of the 
hadron pair is close to the mass of resonances. Examples 
are the p for pion pairs and the for kaon pairs. 

3. A polarimetry 

Another possibility is polarimetry of A baryons 
produced in fragmentation ll2Tll . Here the relevant 
fragmentation function describes the fragmentation 
of a transversely polarized quark into a transversely 
polarized A. The A-polarization can be inferred experi- 
mentally from the kinematics of the decay. 

These three spin-dependent fragmentation processes 
can be used to obtain transversity quark distributions 
1I22II . However, in all three cases the experimental 
observables depend on transversity distributions folded 
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with a currently unknown spin analyzing fragmentation 
function. Therefore transversity functions cannot be 
extracted from polarized semi-inclusive measurements 
alone, and an independent measurement of the fragmen- 
tation functions is required. In this paper we discuss a 
measurement of Collins asymmetries performed at the 
KERB e+e~ collider with the Belle detector. Extraction of 
transversity distributions through a combined analysis of 
transverse spin asymmetries from semi-inclusive polar- 
ized DIS and Collins fragmentation data from e+e^ re- 
quires the universality of the Collins fragmentation func- 
tion between e+e^ and deep inelastic scattering. At 
the present time universality has been confirmed only 
at leading order by Collins and Metz ll23ll . This paper 
describes the measurement of the Collins effect in e+e^ 
collisions. 



quark spin 




FIG. 1 : Definition of the azimuthal angle 4> in the fragmentation 
of a transversely polarized quark with momentum k into an 
unpolarized hadron with transverse momentum P^x. 



B. Collins fragmentation in e+e and SIDIS 

The Collins fragmentation function describes the cre- 
ation of hadrons with transverse momentum P^i^ from 
a transversely polarized quark with polarization and 
momentum k. Following the Trento convention 024(1 the 
number density for finding a hadron h produced from a 
transversely polarized quark q is defined as: 



(k X P^^) ■ S, 

zMh 



(1) 

where the first term is the unpolarized fragmentation 

being the fractional energy of 



CMS 2Eu 

Q 



function, with 

the hadron relative to the center-of-mass system (CMS) 
energy ^/s = Q. The second term contains the Collins 
function H^'^{z, P^j^) and depends on the spin orienta- 
tion of the quark. It changes sign when the quark spin is 
flipped and thus generates a single spin asymmetry. The 
vector product introduces a cos(0) modulation where cf) is 
the azimuthal angle spanned by the transverse momen- 
tum and the plane normal to the quark spin along the 
quark's momentum (see Fig. [l]) . 

In SIDIS one measures a product of the quark trans- 
verse polarization and the Collins function. The trans- 
verse polarization vector of the quarks serves as a refer- 
ence axis for the azimuthal asymmetries of single spin 
distributions. In e+e^ with unpolarized beams no such 
reference direction exists. An analysis of Collins asym- 
metries in single jets in e+e^ therefore will yield a zero 
result as the cosine modulation will average to zero in a 
large event sample. This observation will later be used 
for important tests of the analysis algorithm and possible 
false asymmetries resulting from detector effects (section 
HVl) . The Collins effect can therefore only be observed 
in the combination of two functions each one creating a 
single spin asymmetry. The correlation of quark and an- 
tiquark Collins functions in opposite hemispheres gives a 
product of two cos(0) modulations for the two azimuthal 
angles 0i and 02, resulting in a cos(0i + 02) modulation 



(see Fig. [2]). In e+e annihilation processes these az- 
imuthal angles are defined as 



sgn [n • {(z x n) x (n x P/ii,2)}] 
z X n fix P/ji 2 



X arccos 



|z X n| |fi X Pfti ■ 



(2) 



where z is a unit vector along the z-axis defined by the 
e+e~ beam direction and n is the thrust axis, used as 
an approximation for the quark-antiquark axis (defined 
below) . 



1. Transverse spin components 

In e+e^ annihilation processes with unpolarized 
beams the spin 1 photon has equal contributions from 
the two lepton helicity states | H — ) and | — h). In the 
case a quark-antiquark pair is created with a CMS angle 
of 6* = tt/2 (see Fig. [2]) both lepton helicity combinations 
contribute equally and the transverse polarization of the 
quarks is zero on average. Hence the quark-antiquark 
pair has antiparallel transverse spin components. For fi- 
nite production angles the probability to have antipar- 
allel spins is proportional to sin^ 9. This kinematic de- 
pendence provides a powerful test of the extraction of 
Collins asjnnmetries in e+e^ and will be discussed later. 



C. Azimuthal asjmimetries and cross section 

Two different azimuthal asymmetries in inclusive di- 
hadron production in e+e^ annihilation will become im- 
portant in the course of the analysis. The calculation of 
these will be first described before discussing the cross 
sections. The method already mentioned in the previous 
subsection translates the definition of the Collins func- 
tion (eq. [I]) into the e+e~ — > qq case, where the initial 





FIG. 2: Definition of the azimuthal angles 4>i and (j>2 of the 
two hadrons, between the scattering plane and their transverse 
momenta Phi± around the thrust axis h. The angle 9 is defined 
as the angle between the lepton axis and the thrust axis. 



FIG. 3: Definition of the azimuthal angle <^o formed between 
the planes defined by the lepton momenta and that of one 
hadron and the second hadron's transverse momentum Phi± 
relative to the first hadron. 



momentum of the quark-antiquark pair is known. The 
quark directions are, however, not accessible to a direct 
measurement and are thus approximated by the thrust 
axis. The thrust axis h maximizes the event shape vari- 
able thrust: 



T 



(3) 



where the sum extends over all detected particles. The 
thrust value varies between 0.5 for spherical events and 
1 for tracks aligned with the thrust axis of an event. The 
thrust axis is a good approximation to the original quark- 
antiquark axis as described in Section IIII Al The first 
method of accessing the Collins asymmetry, AI12 is based 
on measuring a cos((/)i + (^2) modulation of hadron pairs 
iN{4>i + 02)) on top of the flat distribution due to the 
unpolarized part of the fragmentation function. The un- 
polarized part is given by the average bin content {N12). 
The normalized distribution is then defined as 



i?12 



N{(t)i -f 02 



{N: 



(4) 



12) 



The corresponding cross section is differential in both az- 
imuthal an gles 4'i>4'2 and fractional energies zi,z2 and 
thus reads 112511 : 



dcr(e^ 



hih2X) 



dVLdzidz2d(f)id(j)2 



'(-l)^?''°'(-2) 



sin2 0cos(0i +02)ifi^'^"'''(2i)i?i''™^''(;^2) ) , (5) 



where the summation runs over all quark flavors acces- 
sible at the center-of-mass energy. Antiquark fragmen- 
tation is denoted by a bar over the corresponding quark 



fragmentation function; the charge-conjugate term has 
been omitted. The fragmentation functions do not ap- 
pear in the cross section directly but as the zeroth ([0]) 
or first ([1]) moments in the absolute value of the corre- 
sponding transverse momenta 1I2F 



Jkrl 

M 



Fiz,]4) 



(6) 



In this equation the transverse hadron momentum 
has been rewritten in terms of the intrinsic transverse 
momentum of the process: P/,j^ = zkr- The mass M is 
usually set to be the mass of the detected hadron, in the 
analysis presented here M will be the pion mass. 

A second way of calculating the azimuthal asjTnme- 
tries, method Mq, integrates over all thrust axis direc- 
tions leaving only one azimuthal angle. This angle is de- 
fined as the angle between the planes spanned by one 
hadron momentum and the lepton momenta, and the 
transverse momentum of the second hadron with respect 
to the first hadron momentum. This angle in the opposite 
jet hemisphere is displayed in Fig. [31 and is calculated as 



ho = sgn [Pk.2 ■ {(z x P,j2) x (P,i2 x P/a)}] 

Z X P,i2 P/i2 X P,,i 



X arccos 



|z X P 



h2\ 



h2 



X Phil 



(7) 



The corresponding normalized distribution i?o, which is 
defined as 



iV(20o) 



(8) 



contains a cos(2(/)o) modulation. The differential cross 
section depends on fractional energies zi , 22 of the two 
hadrons, on the angle 0o and the transverse momentum 
Qt = IqtI of the virtual photon from the e+e^ annihila- 
tion process in the two hadron center-of-mass system. At 
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leading order in as and l/Q^ it assumes the form Il27f1 : 

da{e+e- hih2X) Sa^ 



dQdzidz2(P(lT 



B{y)cos{2cj)o)T 



-zizl\A{y)T{DM 



(2h • kyh • pt — • Pt 



AfiA/2 



>,(9) 



where the convolution integral over the transverse mo- 
menta kr = ziP;ii± and p^ = z2P/i2± is abbreviated 
as: 

T{X\=^el J d^kTd^pTS^{pT + kT-(lT)X . (10) 

a, a 

h denotes a unit vector in the direction of the transverse 
momentum of the first hadron relative to the axis de- 
fined by the second hadron. In Eqn. Q as well as in 
some formulas to follow, the explicit z dependence of the 
fragmentation functions will be omitted while the indices 
will be partially retained: Dl{zi) Di, d\{z2) D2 

and similarly for the Collins functions and H2 ■ The 
kinematic factors are defined as: 



CMS 



{l+COS^e) (11) 



B{v) 



y(i - v) 



CMS 



-{s\T?e) 



(12) 



where = (1 + cos0)/2 is a measure of the forwardness 
of the scattering process. The cross sections given in ([5]) 
and ([9]) are related. Integrating either over the azimuthal 
angles and (j)2 or over the azimuthal angle 0o and 
in the other case will give the same unpolarized cross sec- 
tion. Similarly, the Collins contributions can be related to 
each other, but here, due to the additional convolutions 
of transverse momenta, one either has to know the in- 
trinsic transverse momentum dependence of the Collins 
function or rely on assumptions. A majority of authors 
assume that the Collins function is a Gaussian in kt but of 
different width than the unpolarized fragmentation func- 
tion m, [H, [3^ . Measuring the two types of asym- 
metries one can infer the intrinsic transverse momentum 
dependence of the Collins function. The calculation of 
the asymmetries using a Gaussian kt distribution will be 
discussed further in section IVll 



D. Unlike-sign, like-sign and charged pion pairs 

Favored fragmentation functions describe the frag- 
mentation of a quark of flavor q into a hadron with a 
valence quark of the same flavor. For example, we re- 
fer to the fragmentation processes u — > 7r+ and d 
TT^ (and charge-conjugates) as favored fragmentation 
and to the fragmentation processes u ^ -k" and d 
7r+ (and charge-conjugates) as disfavored fragmenta- 
tion processes. We define yields, iV^ for unlike-sign in- 



clusive pion pair production: 



e ' e 



for like-sign pion production e+e" — > -k'^tt^X and 
N'-^ for charged pions without charge sign identifica- 
tion: e+e^ nnX. Consider for example, the pro- 
duction of unlike-sign charged pions from a pair of up- 
and anti-up-quarks: e+e^ — > uu — > tt^tt^X. The pion 
pair can be either created through two favored fragmen- 
tation processes, D^^ (zi) x (22) or, suppressed in 
yield, through two disfavored fragmentation processes: 
£)J (zi) X D^^{z2). We introduce the favored frag- 
mentation functions d(°-" = Df^^izi) = Dl^ (zi) and 
jjfav ^ x)/ai.(^2) = Dl'{z2) as Well as the disfavored 
fragmentation functions Df'-" = D'^'^'^{zi) ~ [zi) 
and Df-^ = D'^^'^{z2) = (22). Ignoring strange- 
and heavy quark fragmentation into pions and assum- 
ing SU(2) flavor symmetry the cross sections for charged 
pion pair production can be written as: 



^ic^(e^ 



7r±7r=FX) 



dfldzidz2 



-dis 



da{e^ 



dndzidz2 



j:4 



da{e-^ 



D2 



— ( 



COS 



nnX) 



dQ.dzidz2 
2 



cos 



D- 



-fav 



D 



-dis 



(13) 



(14) 



.(15) 



E. Competing sources for azimuthal correlations 

In addition to the Collins fragmentation of back-to- 
back pion pairs there are other processes that also can 
result in azimuthal correlations between the pions. A de- 
tailed understanding of these backgrounds is crucial for 
the extraction of the Collins function. 



1. 'y — Z interference 

Pure Z exchange at Belle energies can be safely ig- 
nored. We briefly consider contributions from 7Z inter- 
ference. Flavor dependent changes to the coupling e^/4 
amount to corrections to the asymmetries of 1.0004 for 
u quarks and 1.001 for d quarks [,27.1 . These corrections 
are small compared to the experimental sensitivity and 
are thus neglected. 



2. Weak decays 

Parity violation in weak decays can lead to azimuthal 
correlations between hadrons in opposite hemispheres. 
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A well known example is the r decay channel r ttv. 
We have used this channel as a test of our extraction pro- 
cedure for azimuthal Collins asymmetries. We have mea- 
sured azimuthal asymmetries in t production and find 
results consistent with existing data llSll. I32l1 and results 
from Monte Carlo. 

In order to remove contributions of t leptons, which 
are characterized by large missing energy, from the 
quark-antiquark data sample we require the energy ob- 
served in the detector to be large, iJvis > 7 GeV. The 
effect of T background remaining after the visible energy 
cut was studied and has been included in the systematic 
error. Similarly, the production of charmed mesons with 
subsequent weak decays can lead to azimuthal correla- 
tions, which have been carefully studied and have been 
subtracted in the extraction of Collins asymmetries (see 
section HVJl ). 



3. Gluon radiative effects 




FIG. 4; Side view of the Belle detector. 



The dominant Collins-like background contribution 
originates from low energetic gluon radiation e+e^ 
qqg, which does not manifest itself in a third jet but cre- 
ates azimuthal asymmetries. The gluon radiation results 
in a cos (2^) modulation and according to calculations 
by Boer 112611 is described by: 



dN 



3 
16^ 



2Q2 



'cos(2(/)o) 



(16) 



The dependence of the asymmetry on is nearly 
quadratic for small and does not depend on the 
charges of the created hadrons. The differential cross 
section in zi and Z2 is proportional to the unpolarized 
fragmentation functions Di{z) and Di{z) for a given 
charge combination of the hadrons. 



II. THE BELLE EXPERIMENT 

The data were taken at the asymmetric e+e^ KERB 
II33I1 storage rings, which collide 8 GeV electrons and 3.5 
GeV positrons. Taking into account the resulting boost 
between the laboratory and CMS frames the Belle de- 
tector itself is asymmetric. It is a large-solid-angle mag- 
netic spectrometer that consists of a multi-layer silicon 
vertex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cherenkov counters 
(ACC), a barrel-like arrangement of time-of- flight scintil- 
lation counters (TOP), and an electromagnetic calorime- 
ter (ECL) comprised of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T mag- 
netic field, see Fig. |4l An iron flux-return located outside 
the coil is instrumented with resistive plate chambers to 
detect mesons and to identify muons (KLM) . The de- 
tector is described in detail elsewhere ll34ll . Two differ- 
ent inner detector configurations were used. For the first 



data sample of 155.6 fb^^, a 2.0 cm radius beam pipe 
and a 3-layer silicon vertex detector were used; the rest 
of the 547 fb^^ of data has been collected using a 1.5 
cm radius beam pipe surrounded by a 4-layer silicon de- 
tector and a small-cell inner drift chamber ifsstl . Particle 
identification in Belle is done by combining data from the 
TOF, ACC and CDC-dE/dx subdetectors. 



A. Monte Carlo 

In order to correct the data for detector effects and for 
systematic studies MC simulated events were generated 
by the Pythia 6.2 or qq98 ll36ll generators and processed 
with a full simulation of the Belle detector based on the 
GEANT3 package lls^ . The MC events are processed 
with the same reconstruction algorithms and written to 
the same data structure as real data. The MC software 
package has a built-in backtracing facility for final state 
particles in order to facilitate the assessment of detec- 
tor effects. In the center-of-mass system the Pythia pro- 
gram generates a back-to-back quark-antiquark pair with 
a 1 + cos^ 9 distribution relative to the e+e^ beam axis. 
The next step is a perturbative parton cascade based on 
order as matrix elements in the leading log approxima- 
tion, which is followed by a parton shower fragmenting 
into the final state hadrons. Within the parton shower 
and hadronization processes transverse momenta are be- 
ing generated, which could lead to azimuthal asjTnme- 
tries. The Collins effect is not modeled in the standard 
event generators used for MC production in Belle. The 
MC samples are divided according to the type of event 
generated: charged and neutral B meson pairs created 
from the T(4S'), light quark (uds) pair production, charm 
quark pair production, and t+t^ production (using the 
Tauola package [38]). In what follows, MC refers to the 
light quark MC simulation if not specified otherwise. 
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III. ANALYSIS 

We report results obtained with an integrated lumi- 
nosity of 547 fh~^. A 55 fb^^ sample was taken at CM 
energy of 10.52 GeV while 492 fb^^ was accumulated on 
the T(4S') resonance at 10.58 GeV. At the lower CMS en- 
ergy, which is below the threshold for BB meson pair 
production, only light and charm quark pair creation 
contribute to the hadronic final states. In the higher en- 
ergy data in addition to continuum events there are res- 
onant T(45) decays into neutral and charged B meson 
pairs. 



A. Event and track selection 

The Collins effect is expected to be dominant in the 
fragmentation of light quarks as helicity is only con- 
served for nearly massless quarks while for heavier 
quarks the correlation between the quark and the anti- 
quark side may be lost. We also focus on the measure- 
ment of the Collins effect in light quark fragmentation, as 
it is the light quark Collins fragmentation function that 
is needed as input for studies of transverse proton spin 
structure in semi-inclusive deep inelastic scattering or 
polarized proton-proton collisions. Most of the B meson 
events can be removed from the data sample, using the 
difference in event shapes between events with underly- 
ing B mesons and light quarks. Since the B mesons decay 
nearly at rest in the CMS, the final state particles exhibit 
a more spherical spatial distribution, which corresponds 
to low thrust values. Most of the light quark-antiquark 
pairs appear in a two-jet topology, which corresponds to 
high thrust values. Consequently, for pion pairs a thrust 
cut of T > 0.8 removes 98% of B data as can be seen in 
Fig. [S] where the simulated thrust distributions for light 
and charmed quark pairs and T{4S) decays are shown. 

For the calculation of the thrust variable all charged 
tracks and all neutral particles with a minimum energy 
of 0.1 GeV are considered. For the purpose of obtain- 
ing an unbiased data sample one assigns the sign of the 
thrust axis at random. The contribution from B mesons 
to the observed asymmetries can be estimated by com- 
paring the data taken on the T(4S') resonance with the 
data taken 60 MeV below the resonance. This test will 
be discussed in Section HVNI Events with charm quarks 
do not exhibit a very different event shape from light 
quark events, see Fig.[5j However, the contributions from 
events with charm quarks can be corrected by measuring 
azimuthal asymmetries in a charm-enhanced data sam- 
ple. This will be described in Section HVJl 

In order to ensure a two-jet geometry in the selected 
event sample with the majority of final state particles re- 
constructed in one of the two jets a minimum visible en- 
ergy of Eyisibie > 7 GeV is required. Charged tracks 
used in the analysis are required to originate from the 
interaction point and to lie in a fiducial region —0.6 < 
cos(0iab) < 0.9, where 6\s.h is the polar angle in the lab- 
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FIG. 5: Simulated thrust distributions for selected 2-pion pairs 
at ^/s = 10.58 GeV, for e+e^ ^ events (open dia- 

monds), e'^e~ —* B'^B° events (open circles), e+e^ cc 
events (full triangles) and for light quark production e+e^ 
qq, q G uds (full squares) normalized to the total number of 
events in all channels. The vertical line represents the minimal 
thrust value selected for the analysis. 



oratory frame relative to the direction opposite the the 
incoming positron (definition of the z-axis). This cor- 
responds to a nearly symmetric fiducial interval in the 
CMS frame -0.79 < cos6'cms < 0.74 and covers the ac- 
ceptance of the barrel part of the Belle detector. For the 
identification of pions a likelihood ratio is used, which 
is based on energy loss in the drift chamber (CDC), the 
number of Cherekov photons (ACC) and time of flight 
information (TOF). Kaons are separated from pions by 
requiring C{tt)/[C{K) + C{t:)] > 0.7. C{tt/K) is the like- 
lihood for a track be a pion or kaon. The percentage of 
misidentified pion pairs is below 10% in all zi and Z2 
bins. In addition, the likelihood ratios for being either a 
muon or an electron have to be below 0.9 and 0.8, re- 
spectively. A cut on the fractional hadron energy of the 
two hadrons zi_2 = 2Ei-2/Q > 0.2, avoids contributions 
from decays with the decay products incorrectly recon- 
structed in opposite hemispheres. 

The two pion tracks are required to lie in opposite 
hemispheres with the selection VFiiomi (P/ii • n)(P;j2 • 
n) < 0, where the hemispheres are separated by the 
plane normal to the thrust axis n. A comparison of 
the quark-antiquark axis with the thrust axis calcula- 
tions from reconstructed particles shows an average an- 
gular deviation between the two of 128 ± 82 mrad (the 
RMS value is quoted for the uncertainty) in simulated 
events for light quark production, while it appears to 
be slightly larger for charm production (see Fig. [6l) due 
to semileptonic decays. Since the thrust axis calculated 
from generated particles also deviates from the original 
quark-antiquark axis by a similar magnitude, we con- 
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FIG. 6: Angular deviation from the original quark-antiquark 
axis of the reconstructed thrust axis (solid line) and of the 
thrust axis calculated from generated events before detector 
simulation (dashed line) . The distributions are based on a light 
quark MC simulation with a T > 0.8 requirement for the gen- 
erated events. 
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FIG. 7: Hemisphere variable W^hcmi compared to the product 
of energies of the two hadrons (top) and the transverse vir- 
tual photon momentum Qt (bottom plot) for pion pairs with 
Whcmi > —5 GeV^. The horizontal line corresponds to the 
maximum Qt value selected for the analysis. 



elude that the observed deviation arises from the intrin- 
sic difference between the original quark direction and 
thrust axis in the hadronization process. A test with the 
approximate algorithm buildt in Pj^hia confirms these 
deviations. The smearing of the thrust axis compared 
to the quark axis leads to a reduction in the amphtudes 
of the measured azimuthal back-to-back correlations, as 
discussed below 

The impact of initial state radiation (ISR) from the 
incoming leptons was estimated using MC simulation. 
The average CMS energy for events fulfilling all selec- 
tion criteria is reduced by ISR to 10.51 GeV for the on- 
resonance data with about 2% of events at energies be- 
low 9 GeV. Hence the fractional energies are slightly un- 
derestimated, but the effect of ISR is negligible compared 
to the width of the z-bins used in this analysis. The 
effect of ISR on the asymmetry measurements was in- 
ferred from a MC simulation employing a reweighting 
technique, which is discussed in section IIVDI The re- 
sults were consistent with those obtained without ISR 
and thus, the overall effect of QED initial state radiation 
on this analysis is negligible. 

Possible hemisphere misassignment is further reduced 
by the requirement that the transverse momentum of the 
virtual photon in the two-pion CMS (Qt) be smaller than 
3.5 GeV. The high Qt region is mainly populated by pion 
pairs with very asymmetric fractional energies z and es- 
pecially when the lower energetic particle happens to be 
close to the hemisphere boundary (see Fig. [T]) . For such 
pairs the remaining uncertainty in the determination of 
the thrust axis direction can lead to a hemisphere misas- 
signment and the simple correlation of one pion originat- 



ing from the quark side and the other pion originating 
from the antiquark side will be spoiled. At the hemi- 
sphere boundary asymmetric particle decays could also 
lead to the assignment of one pion into the wrong hemi- 
sphere. After application of the Qx cut the hemisphere 
is misassigned in less than 0.1% of the events according 
to MC simulation. 



B. Normalized jdelds, raw asjonmetries 

In the expression for the di-hadron yields in ^ and 
(O the product of two Collins functions appears as the 
amplitude of a cosine modulation in the di-hadron yield. 
The Collins functions depend on the fractional energies 
zi.2 of the two hadrons. Therefore the analysis is per- 
formed in 4 bins in zi and Z2 for each hadron with 
boundaries at Zi = 0.2, 0.3, 0.5, 0.7 and 1.0. Alterna- 
tively, the cosine modulation will be evaluated in bins of 
Qt- In this case the dependence on the fractional en- 
ergies is averaged over all bins in zi and 2- As a first 
step in the analysis, the two-pion yields NaiPa) are ob- 
tained for each (2:1,22) bin in 8 equidistant bins of the 
azimuthal angles a = 0, 12. The index a = 12 refers 
to the thrust method M12 introduced in ([5]) and the angle 
/3i2 ^ (1)1 + 4'2 is the sum of the hadron angles shown in 
Fig. [21 The index a ~ refers to the two hadron method 
Ma as introduced in Eqn. ^ and the angle po ~ 20o is 
defined in Fig. [3j Next, the azimuthal di-hadron yields, 
NaiPa) are normalized to the average di-hadron yield 
(Na) per bin resulting in the normalized yields i?^, de- 
fined in Eqns. ^ and dS]). Following Eqns. Q and dH) 
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FIG. 8: Raw data. Top: Example of uncorrected unlike-sign 
(open circles) and like-sign (open squares) pion pair normal- 
ized rate Rq vs. 2(f>o in the bin 21(22) G [0.5,0.7], 22(21) G 
[0.3, 0.5]. Bottom: Pion pair double ratio Rq /Rq vs. 2(^o in the 
same 21 , 22 bin. The results of the fit described in the text (full 
and dashed lines) are also shown. 
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FIG. 9: Raw asjnnmetry parameters for unlike-sign pion 
pairs and for like-sign pairs as a function of the transverse 
virtual photon momentum Qt- Light quark MC generated (tri- 
angles), reconstructed (squares), and selected data sample (cir- 
cles) for the Q = 0, 12 reconstruction methods are shown. 



asymmetries for pion pairs from generated and recon- 
structed MC data samples. All track and event selection 
cuts have been applied with the exception of the cut on 
the transverse photon momentum Qt < 3.5 GeV. No 
particle identification and polar angle selection cuts have 
been applied for the generated MC sample. 



C. Double ratios 

The flo raw asymmetries obtained from generated MC 
events are nonzero and increasing with Qt while the ai2 
are almost zero. A large difference between the asjnnme- 
tries calculated directly from generated MC events and 
the asymmetries obtained from MC events after GEANT 
simulation of the detector response and track reconstruc- 
tion can be observed in Fig. [9) This difference points to 
large acceptance effects, which can even cause the Aq 
asymmetries to change sign in method AIq. We have 
studied the possibility of cancelling acceptance effects in 
the asymmetries oq and ai2 using a MC simulation of the 
acceptance effects but found it difficult to separate asym- 
metries from the Collins effect from acceptance effects 
and the radiative background in Oq. The small differ- 
ences between the asymmetries from reconstructed 
MC events and real data gives an estimate on how well 
radiative effects are described in the MC event genera- 
tor Therefore in order to correct the raw asymmetries 
for detector effects one has to rely heavily on a MC simu- 
lation, which does not include any spin-dependent asym- 
metries. As a consequence we do not consider these any 
further and employ a scheme using ratios of asjnnme- 
tries, in which most of the instrumental effects cancel, 
and is thus independent of the details of the MC simula- 
tion. 



1. Method 



the normalized yields can be parameterized as: 

Ra = cos(/3q) + 6q, (a = 0, 12). (17) 

The measured normalized yields Ra are fitted with the 
parameterization in ( |17t leaving Cq and ba as free pa- 
rameters. The constant terms ba are found to be consis- 
tent with unity within the statistical errors as shown in 
Fig-in The raw asymmetry parameters, Oa, are found to 
be large as can be seen in Fig. [9] for the unlike-sign (U) 
and like-sign (L) pion pairs as a function of Qt- The fig- 
ure compares asymmetries for pion pairs from data with 



The asymmetries generated by QCD radiative events 
and acceptance effects do not depend on the charge com- 
bination of the pion pairs. For detector effects this can be 
tested by comparing positively charged with negatively 
charged pion pairs, as will be discussed in section IIVFI 
When building the normalized yields Ra the normaliza- 
tion is obtained by integrating over the azimuthal angles 
(3a- This causes the terms proportional to cos(/3q) to van- 
ish. Adding to the cross section in (O the term due to 
gluon radiation, which is proportional to the unpolarized 
fragmentation functions, and has been scaled by a factor 
C, and dividing by the average number of hadron pairs 
in the full Pa range, we obtain expression (I18D for R12- 
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(18) 



In order to obtain a MC independent measure of the 
Collins effect, one can form a double ratio of, for exam- 
ple, the normalized yields of unlike-sign over like-sign 
pion pairs: 



Nkm/{Nk) 



0,12) 



(19) 



Up to the linear term in the amplitude of cos(/3q) the 
term multiplied by C in ( flSl) cancels and one is left with 
the following expression for the double ratio: 
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(20) 



Similarly one obtains the double ratio of unlike-sign and pairs: 
any sign (i.e. ++, H — and charge-conjugates; C) pion 
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This double ratio contains a different combination of fa- 
vored and disfavored fragmentation functions as pointed 
out by Efremov and Schweitzer 113011 . Analogous expres- 
sions can be obtained for the Rq/Rq double ratios with 
i,j € U,L, C. In the double ratios the acceptance ef- 
fects cancel while the QCD radiative effects cancel to 
first order. The possible influence of higher order terms 
(cos^(/?q)) has been studied explicitly by including these 
terms in the fits and also by performing a subtraction of 
U and L (C) normalized yields, as discussed in the next 



section. Statistical correlations between the U, L and C 
pairs have been taken into account. The double ratios 
are parameterized by R^/Ra — cos(/3ct) + B^J, and 
the measured distributions are fitted with A^^ and B^J as 
free parameters. The constant terms B^^ are again found 
to be consistent with unity within statistical errors. 
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IV. SYSTEMATIC STUDIES 



B. Single spin asymmetries 



The impact of the detector performance and of the 
method used for asymmetry reconstruction on the iso- 
lation of the CoUins effect has been estimated through 
various systematic studies. 



A. MC double ratios 

An important test of the double ratio method is the 
extraction of double ratios from MC. The generic MC de- 
scribes the radiative gluon effects. The acceptance effects 
are also included in the detector simulation. However, 
the MC generator does not contain asjnnmetries based 
on the Collins effect. Therefore one expects the cosine 
moments of the double ratios for MC to vanish. Fig- 
ure [To] shows the fitted asymmetry parameters as 
a function of the zi, Z2 bins. It can be seen that these 
asymmetries are consistent with zero in all bins. Similar 
features are observed for the asymmetry A^^ , which is 
not shown. A fit of a constant to the observed asymme- 
tries Aj^^ and A^"^ was performed and the largest devia- 
tions from zero, including the statistical errors of the fits, 
were attributed as systematic uncertainties of the final 
results. The contribution to the systematic uncertainty 
ranges from 0.07% in to 0.22% in A^'^ as summa- 
rized in Table [nil 
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Aq^ (circles) and 

j4 (squares) as a function of 22 for 4 21 bins obtained 
from MC double ratios. The results of fits to a constant are 
also displayed as dash-dotted lines. They are (0.13 ± 0.09)% 
for A'i^, (0.08 ± 0.09)% for A'd' , (0.06 ± 0.04)% for AY'^ and 
(0.04 ± 0.04)% for AY2 ■ The UC data points are not shown. 



The Collins fragmentation in e+e^ annihilation ap- 
plies when the two primary quarks have equal but oppo- 
site transverse spin components; asymmetries only occur 
when hadrons from both quarks are observed simulta- 
neously. For unpolarized beams there is no net trans- 
verse polarization of an individual quark or antiquark, 
and hence there is no modulation of 0i or 02 separately. 
Using the same event selection, and measuring sin 4>i or 
sin 02 modulations we found these asymmetries to be 
consistent with zero within the statistical uncertainty. No 
systematic error is assumed. 



C. Event mixed asjTiimetries 

The spin information of quark-antiquark pairs also is 
lost if one uses the same analysis procedure but com- 
bines two pions from different events. Such pion pairs 
would be uncorrelated and hence no asymmetries should 
be visible. This test has been performed and the results 
were consistent with zero asymmetries with the largest 
asymmetries being below 0.01%. No contribution to the 
systematic error has been assigned. 



D. Moment reconstruction in reweighted simulated 
samples 

The MC generators do not include spin effects, there- 
fore in order to perform a realistic simulation, which 
is as close as possible to the experimental observations, 
weights have been introduced to generate asymmetries. 
This allows to study the influence of the detector effects 
on the asymmetries. The weights are applied to the 
generated particles before detector simulation and de- 
pend on the azimuthal angles /3q and the fractional en- 
ergies Zgen of the final state particles. 



= 1 + a\Zgen) cos (3a,gen , {i = U, L, C; 



0,12) 



(22) 

Analyzing the UL double ratios with, for example, gen- 
erated weights = 5% and = —5%, the recon- 
structed double ratios should ideally return a 10% asym- 
metry. An example of the reconstructed UL asymme- 
tries as a function of the z variables is shown in Fig. 1111 
compared to the linear weights in zi and Z2 of = 
1 + 0.05 • Z1Z2 cos(/3q) and ~ \ — 0.05 • ziZ2 cos(/3q,). 
A summary of several combinations of zi 2 independent 
weights and their reconstructed values (as obtained by a 
constant fit to all z bins) are displayed in Table HI The 
Aq results are consistent with the generated asymme- 
tries, while the A12 results systematically underestimate 
the generated asymmetries. This underestimation can be 
attributed to the difference between the reconstructed 
thrust axis and the original quark-antiquark axis. The di- 
rection of this axis is essential in the estimation of A12; 
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FIG. 11: Double ratio results Aq^ (circles) and AYi^ (squares) 
as a function of Z2 for 4 zi bins for reconstructed light quark 
MC, reweighted as described in the text. The generated asym- 
metry Q.lziz-z evaluated at the bin centers is also shown as his- 
togram. The dotted line corresponds to the zero asymmetry. 



because of the deviation the calculated azimuthal angles 
and also their modulations are smeared and thus the am- 
plitudes of the modulations decrease. An average un- 
derestimation of (58.9 ± 1.3) % in the A12 results is ob- 
served. In order to obtain asymmetries relative to the 
quark-antiquark axis these results are thus rescaled by a 
factor 1.66 ± 0.04. The corrected asymmetries can now 
be used at all energies, because the deviation between 
the thrust axis and the quark directions, which might 
be energy dependent, is removed. The Aq asymmetries 
depend only marginally on the original quark-antiquark 
axis due to the hemisphere assignment. The small un- 
derestimation of the asymmetries can thus be attributed 
to the smearing of the tracks in the detector. To correct 
for this dilution, the Aq results are rescaled by a factor of 
1.11 ± 0.05. The error on the correction factor is added 
as systematic error (see Table HiH). 



TABLE I: z-independent weights generated in the U, L and 
C channels using generated tracks and the original quark- 
antiquark axis and the reconstructed average asymmetries as 
obtained by a constant fit to all z bins. 

weights (%) UL double ratios (%) UC double ratios (%) 

au aL ac Aq A12 Aq A12 



+5 -5 
10 



-5 




9.8 ±0.1 
9.8 ±0.1 



5.9 ±0.2 
5.8 ±0.2 



9.8 ±0.1 
9.8 ±0.1 



6.0 ±0.2 
5.9 ±0.2 



E. Double ratios and subtraction of normalized yields 

The method of building double ratios, as described 
above, cancels possible acceptance effects as well as ra- 



diative effects to leading order However, higher orders 
in the expansion of the radiative term might still remain 
and could affect the results if they were large. A second 
method exists that will cancel the radiative terms, but 
not necessarily the acceptance effects. If one subtracts 
the normalized yields for one charge combination from 
those of another charge combination. 
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one is sensitive only to the Collins asjmimetry and pos- 
sible acceptance effects. Again these yields are fitted by 
a constant and a cosine modulation 5'^^ = A^^ cos(/3a) ± 
B^J, with i,j € U, L, C where the constant fit parameter 
now is consistent with zero. As a systematic study 
the results obtained from the double ratio method and 
the subtraction method have been compared and the ab- 
solute value of their differences has been assigned as a 
systematic error on the reconstruction method. As sum- 
marized in Tablellnl these differences are small. The mag- 
nitude of the systematic uncertainty in the Aq^ results 
ranges from less than 0.01% in the bins with small frac- 
tional energies to 0.26% in the highest z bin where the 
smaller statistics might drive the differences. 



F. tt+tt+Ztt tt double ratio tests 

Another possible source of systematic errors could be 
a charge dependence of the detector response, which 
would manifest itself in azimuthal asymmetries. One can 
test this by probing double ratios of positively charged 
pion pairs over negatively charged pion pairs. Neither 
the Collins asymmetry nor the radiative effects depend 
on the pion charges. Therefore, one expects zero asym- 
metries. Again we fit the Aa results for all z bins by a 
constant and assign the result together with its statistical 
error as a systematic uncertainty. The magnitude of the 
absolute systematic uncertainty is 0.05% as summarized 
in Table [IlTl 



G. r+r contribution to the asymmetries 

Background from r pairs, which survive the selection, 
may give rise to fake asjnnmetries. To estimate this effect 
one has to use a data set that contains mostly e+e^ 
T^T^ events. The t content is enhanced by requiring a 
visible energy less than 10 GeV. In this sample about one 
third of the events originates from r pairs while the tau 
contributions in the real data for the combined z binning 
(symmetric zi,z2 were combined into one bin) are: 



bin 



1 



4 5 6 7 8 9 



#7r^/#7ra;, 1% 2% 2% 4% 2% 3% 5% 4% 7% 8% 
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When analyzing the r enhanced data in the same way 
as the full data sample one finds the following asjTnme- 
tries: A^^ = (0.179 ± 0.300)% for the cos2(^o method 
and A^^ = (0.750±0.300)% for the cos((/)i +^2) method. 
These values are consistent with zero within 2a but con- 
sistently slightly positive; these results can be explained 
by the fairly large residual e+e~ ^ qq,q € udsc contri- 
bution in the t enhanced data sample. As a consequence 
it is not necessary to correct the measured asymmetries 
for this contribution. No systematic uncertainty has been 
assumed. 



H. Uncertainties due to particle identification (PID) 

While the applied pion selection only admits less than 
10% of misidentified pion pairs, the possible contri- 
butions from this misidentified background has to be 
checked. For this purpose a tighter pion likelihood selec- 
tion has been applied, where the pion-kaon separation 
requirement was tightened to £(7r)/[£(_ft') +£(7r)] > 0.9. 
No significant changes have been observed. The small 
differences were added to the systematic uncertainties, 
which range from 0.01% to 0.05% (see Table |llll). 

I. Uncertainties due to higher harmonics in the fit 

Furthermore one has to ensure the robustness of the 
fit. For this purpose one compares the cosine moments 
taken from the fit described above with the cosine mo- 
ments when including an additional sine modulation and 
a cosine modulation of twice the argument as free pa- 
rameters. No significant changes are observed and there- 
fore no systematic uncertainty has been assumed. In ad- 
dition, varying the number of bins in azimuthal angle pa 
between 6 and 24 gave stable results for the extracted 
cos f3a moments. 

J. Charm contributions and correction 

The contribution from e+e~ — *■ cc amounts to about 
40% of the total quark-antiquark production cross- 
section. Due to weak decays of charmed hadrons, which 
can introduce azimuthal asymmetries not originating 
from the Collins effect, the results have to be corrected 
for this contribution. 



1 . Charm enhanced data sample 

The decays of charmed hadrons are well described 
in the MC and, one can derive the relative contribu- 
tions D := [data of pion pairs in the data sam- 
ple. Additionally one can select a charm-enhanced data 
sample, where one requires a D-k pair invariant mass 
(m{D'K)) consistent with the mass of the D* meson. 



corresponding to decays D*~^ D°7r+ (and charge- 
conjugated process). These decays are selected with a 
high purity by requiring (in addition to the selection 
criteria already described) the invariant mass difference 
Am = m{Dn) — m{D) consistent with the nominal mass 
difference of D* and D mesons its^ . One can calcu- 
late the charm quark-antiquark pair contributions to this 
data sample as d (see Table [III). Measuring 

the double ratio asymmetries A^^ simultaneously in the 
default as well as in the D* meson enhanced data sam- 
ples and assuming that the asymmetries from the charm 
events are the same in both samples, one obtains two 
equations for the two separate data samples: 

^uds = ~^ J^Agii — ^^charm (25) 

^charm = ~,^D' J — ^uds j (26) 

d d 

where ^aii are the measured asymmetries in the default 
data sample. Ad* those of the charm enhanced data sam- 
ple and Auds and ^charm the true asymmetries of uds 
or charm quarks, respectively. As a result one obtains 
corrected asjnnmetries for the light quarks. The statis- 
tical errors of both samples were propagated to obtain 
the statistical errors for A^^s (and Achaim)- This leads to 
comparatively large statistical errors for the charm cor- 
rected asymmetries due to the low statistics in the charm 
enhanced data sample. 

K. Bottom background 

Due to the thrust selection and the requirement on the 
fractional energy of the tracks only a small fraction of 
pion pairs from BB decays appears in the data sample. 
Only the bins with fractional energies below 0.5 can have 
a B -contribution. As a consequence the contributions 
B := ^if^ is at most 2.5% as can be seen from Table 

iVall 

ini Correspondingly the systematic uncertainties due to 
B meson decays in data are negligible compared to other 
uncertainties. In correcting the light quark asymmetries 
for the charm contribution the small contributions from 
bottom background were taken into account by replacing 
1 — Dhyl — D — B and similarly 1 — dhyl — d — b. The 
bottom quark asymmetry has been assumed to be zero. 

L. Beam polarization measurements 

Another important study is a test of the polarizations 
of the electron and positron beams. The natural beam 
polarization in electron storage rings is a transverse po- 
larization with the electron spins aligned with the mag- 
netic bending fields. This is a consequence of the emis- 
sion of synchrotron radiation as detailed in a paper by 
Sokolov and Ternov ll40tl . In the KERB storage ring the 
polarization will be destroyed by the large tune shifts 
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TABLE II: Relative charm, B meson and uds contributions for the 2:1,2:2 binning in the selected data sample and the D* sample in 
%. 
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in the beam-beam interactions. A measurement of any 
residual polarization is, however, in order. While the pro- 
cess e+e^ 77 appears to be the most sensitive to a 
transverse beam polarization, the process e+e^ M^A* 
is experimentally easier to asses. For muon pair events 
the cos{d) and (t>LAB distributions are analyzed, where 
9 is the polar angle between the electron and the muon 
axis in the CMS and (pLAB is the azimuthal angle of the 
muon in the laboratory frame around the z-axis. Both 
distributions are normalized to the average bin content. 
In the case of the cos{6) distribution this procedure was 
performed only in the range [—0.75, 0.75] since at larger 
angles acceptance effects dominate. The cos(0) distribu- 
tion is then fitted with the functional form a(l + 6 cos^ 9 + 
ccos6). The unpolarized cross section should contain a 
cos^{9) term equal to unity and a cos{6) moment of a 
few percent, due to photon-Z interference and higher or- 
der QED terms. An example of these fits is shown in 
Fig. [12I The 6 moments show the expected modulations 
of = 1.02 ± 0.04 and c = 0.056 ± 0.0015 and show that 
the muon selection and detector performance is reliable 
for muon pairs. The (I)lab distribution is fitted by the 
function a + 6cos((?!)lab) + csin(^LAB) ± dcos{2(j)LAB)- 
An example of such a fit is shown in Fig.ll2l Any nonzero 
(j)LAB moment would be a sign of a nonzero beam polar- 
ization. The cos2cj)LAB moment is especially sensitive 
since it is proportional to the product of the two beam 
polarizations Ii41]. The fits to the 20lab distributions 
were performed either by taking all the muon pairs of 
one run or by dividing individual runs into several time 
bins to study the time behavior of a possible beam po- 
larization. Run periods before and after the implemen- 
tation of continuous injection as a mode of KEKB opera- 
tion Ii33i1 were also compared. The c/jlab moments are 



consistent with zero in all time periods. No build-up of 
polarization could be observed. Therefore one can con- 
clude that no significant beam polarization exists in this 
data sample and no systematic error is assigned. 

M. Correlation studies 

The statistical errors obtained from the fits to the dou- 
ble ratios may be influenced by correlations among indi- 
vidual angular and (zi, Z2) bins. The correlations can 
arise since different pairs composed of the same pion 
are allowed to contribute to the distributions. The effect 
on the statistical errors was tested by dividing the large 
reweighted MC sample into a large number of subsam- 
ples. The same fit as for the data was performed using 
each of the MC subsamples. The width of the distribu- 
tion of the fit results is found to be in agreement with 
the statistical error of the fit performed on the full MC 
sample. 

N. Internal consistency test 

Before combining the published data f3] obtained 
from the off-resonance sample and the additional on- 
resonance data, the consistency of the two data sets 
has to be tested. The charm corrected results from 
the published 29.1 fb^^ data sample and the results 
from the 492 fb~^ data are displayed in Fig. 1131 the 
results are in good agreement in all combined z-bins. 
For this comparison the correction factors for the pub- 
lished data, which are based on the generated thrust 
axis instead of the quark axis were applied. The overall 
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FIG. 12: Angular distributions of cos^cms (top plot) and 
2cj>LAB (bottom plot) for muon pairs. 



X 



{Ai{setl)-Ai{set2)y 



of the two data sam- 



1,10 AA'f(setl)+AA'f{set2) 

pies per degree of freedom is 0.54 for the Aq^ results 
and 0.83 for the Afj^. Therefore it is possible to combine 
the off-resonance and on-resonance data sets. The differ- 
ences in the CMS energy are automatically accounted for 
by considering fractional energies z. 

An additional test to compare the results of differ- 
ent data taking periods for on-resonance data and off- 
resonance data and verify that these are consistent as 
well as independent of time. At Belle the data is natu- 
rally divided into periods of several-months' data taking, 
called experiments with odd, increasing numbers start- 
ing at 7. As a reference the double ratio results for the 
complete data set (experiments 7-49) for both types of 
double ratios ( UL and UC) and extraction methods (0 
and 12) have been taken. For each experiment number 
the value per degree of freedom relative to the com- 
bined result is calculated. The values are presented 
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FIG. 13: Comparison of the A'^^ asymmetries for the published 
off-resonance data sample (triangles) and the on-resonance 
data sample (squares). The upper plot shows the Ao asym- 
metries, the lower plot the A12 asymmetries as a function of 
z. 



in Fig. [14] as a function of the experiment number No 
systematic trend for any of the data samples or methods 
can be observed. The distributions of x^ per degree of 
freedom are also displayed in Fig. I15[ where one sees 
that they scatter around the central value of 1. 
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FIG. 14: distribution per degree of freedom of the A^^ (tri- 
angles) and A^'"' (circles) asymmetries as functions of the ex- 
periment number for the Aq (top) and the A12 (bottom) asym- 
metries. The open symbols represent the off-resonance data 
samples, the full symbols the on-resonance data samples. 



In summary a number of possible sources of uncertain- 
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FIG. 15: distribution per degree of freedom of the A'^^ (top 
plots) and A^"^ (bottom plots) asymmetries for the Ao (left 
plots) and the A12 method (right plots). 



ties in the asjnnmetry extraction have been studied and 
their contributions have been evaluated. As can be seen 
in Fig. [16] and Table [Hi] the errors are dominated by the 
detector effects on the double ratios and the statistical 
uncertainties on them. In general, most of the systematic 
uncertainties have significantly decreased in comparison 
to the previously published data Ist] as the statistics of 
the data that are used to evaluate some of the systematic 
uncertainties also increased by a factor of almost 20. 



V. RESULTS 

The final results combine the 55 fb~^ data sample 
taken at an energy of 10.52 GeV and 492 fb ~^ of data 
taken on the T(4S') resonance at 10.58 GeV. Since the 
fractional energies zi,2 are already normalized by the 
corresponding CM energies the two data sets have been 
combined. The double ratios have been evaluated and 
a fit was performed as described above. The asjnnme- 
tries have been corrected for the charm contribution and 
were rescaled by the factors obtained by the weighted 
MC (1.66 ± 0.04 for the asymmetries, 1.11 ± 0.05 for 
the Ao asymmetries) . 



A. Double ratio results 

1. Double ratios versus fractional energies z\zi 

The main results are the asymmetry parameters ^0 
and A\2 for both types of double ratios (UL and UC) as 
a function of the fractional energies of the two hadrons. 
Figures [17] and [18] show these asymmetries where all 
bins for a given z\ are displayed. The numerical values 
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FIG. 16: Squares of the contributions to the systematic error 
for the A^-t asymmetries as a function of z\ and Z2. The lower- 
most bars correspond to the contribution from the differences 
between the double ratio and subtraction reconstruction meth- 
ods. The next bars correspond to the charged pion pair ratios, 
the next to the MC double ratios, then the contribution due 
to particle identification and the uppermost bar corresponds to 
the uncertainties in correcting the measured asymmetries for 
the quark-antiquark/ thrust axis deviation. 



are give in Tables [W] and (V] One can clearly see the rising 
asymmetry in each plot as a function of Z2- The UC asym- 
metries are significantly smaller than the UL asymmetries 
but non-zero, which, given the different contributions of 
favored and disfavored fragmentation functions already 
suggests a large disfavored Collins fragmentation func- 
tion with opposite sign to the favored one. This sugges- 
tion will be quantified in the next section. 



2. Double ratios versus polar angle sin^ 6/(1 + cos^ 9) 

Another way of presenting the results of the asjnnme- 
try measurements is based on the requirement for pri- 
mordial transverse quark polarization. At leading order 
and ignoring 7 — Z interference transverse quark spins 

leads to a ^^^'"pgf^g^ dependence of the asymmetries (see 
Eqs. 1201 1211) . where 6 can be represented by either the 
polar angle between the e+e^ and the thrust axis n^, or 
by the polar angle 62 of the 2"'' hadron relative to the 
e+e~ axis. Figure [19] displays the Aq and A12 results 



for the UL double ratios as a function of , .> while 
Fig. [20] displays the results for the UC double ratios. Both 
polar angles have been considered and each of them has 
been fit by a first order polynomial where the constant 
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TABLE III: Absolute contributions to the sj^tematic error on the and A^ results for the zi, zi binning. The different columns 
contain the systematic errors arising from subtraction method (1), ratios of positively over negatively charged pion pairs (2), 
double ratios from MC (3), uncertainties due to particle identification (4) and due to underestimation in reweighted MC (5) as 
described in the text. 
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FIG. 17: Light quark (uds) Ao asymmetry parameters as a func- 
tion of Z2 for 4 2i bins. The UL data are represented by triangles 
and the systematic error by the upper error band. The UC data 
are described by the squares and their systematic uncertainty 
by the lower error band. 
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FIG. 18: Light quark (uds) An asymmetry parameters as a 
function of 22 for 4 zi bins. The UL data are represented by 
triangles and the systematic error by the upper error band. The 
UC data are described by the squares and their systematic un- 
certainty by the lower error band. 



term has been set to zero. In both cases the results are 
consistent with a linear behavior The results obtained 
with the thrust axis defining the polar angle can be de- 
scribed by the linear term only as the per degree of 
freedom of the fit changes only slightly when allowing 
the constant term to float, for example for the A^^ re- 
sult from 2.4 to 1.67 and from 2.56 to 2.35 for the 
result. The Aq results obtained with 62 as the polar angle 
favor a nonzero constant term; when a constant term is 
included the reduced of the fit decreases significantly 
from 2.81 to 1.26 for the A^'^ resuh and from 2.57 to 
1.22 for the A\^^ result. This can be explained by the 
fact that the thrust axis describes the original quark di- 
rection better than the 2"'* hadron's polar angle, which 
receives some additional transverse momentum relative 
to the quark axis. 



3. Double ratios versus Qt for high and low thrust data 
samples 

The dependence of the asymmetries on the virtual pho- 
ton momentum in the two-hadron center-of-mass frame 
is also of interest. The results are shown in Figs. [21] 
and 1221 In addition to the charm-corrected asymmetries 
the asymmetries for the reverse thrust selection T < 0.8 
are displayed. The contributions of both charm quarks 
and by T (45*) decays are quite substantial in the reverse 
thrust selection sample and can add up to almost 70% 
in the highest Qt bin. The results of the reverse thrust 
selection are displayed uncorrected for the charm and 
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FIG. 19: Light quark (uds) (top) and (bottom) asym- 
metry parameters as a function of sin^ S/(l + cos^ ff), for Q2 
(squares) and for fiz (triangles). Linear fits are also displayed 
as dashed and continuous lines, respectively. The systematic 
error for Q2 case is represented by the lower, that for by the 
upper error band. 



the T(4S') contributions. When comparing the reverse 
thrust selection for on and off-resonance data one sees 
that the T(45') does give an additional contribution to 
the A12 asymmetries. Nevertheless it is clearly visible 
that the asjnnmetries are significantly lower than in the 
main data selection. This is the expected behavior, since 
the asymmetries due to the Collins effect are smeared out 
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FIG. 20: Light quark (uds) Aq (top) and A^f (bottom) asym- 
metry parameters as a function of sin^6/(l + cos^6'), for 62 
(squares) and for (triangles). Linear fit are also displayed 
as dashed and continuous lines, respectively. The systematic 
error for 62 case is represented by the lower, that for by the 
upper error band. 
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rected for heavy quark contributions (squares). The vertical 
line represents the main data selection Qt < 3.5 GeV. 
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FIG. 21: Light quark (uds) Aq^ (top) and A^'a^ (bottom) asym- 
metry parameters as a function of Qt, for events with T > 0.8 
(triangles) and asymmetries for events with T < 0.8 not cor- 
rected for heavy quark contributions (squares). The vertical 
line represents the main data selection Qt < 3.5 GeV. 



for events with no clear two-jet structure. 



4. Charm asymmetries 

Eqns. (I25D and (1261) can also be solved for the asym- 
metries from the process e+e^ cc. The results as a 
function of the fractional energies zi and Z2 are displayed 
in Fig. [23] While a small asymmetry is visible at lower 
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FIG. 23: Charm asymmetry parameters ^0 and A12 as a func- 
tion of the combined z. The UL data are described by triangles, 
their systematic error being the top error band while the UC 
data are described by the squares and their systematics by the 
lower error band. 



fractional energies it seems to be consistent with zero at 
larger fractional energies, although the statistical errors 
become large. The integrated results for the Aq asymme- 
tries are compatible with zero i{A^'^) = -0.011 ± 0.007 
and {A^'^') = -0.004±0.003), while the A12 asymmetries 
are found to be slightly positive ({AY2) = 0.037 ± 0.011 
and (A^^) ^ 0.012 ± 0.003). 
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VI. PARAMETERIZATION OF THE COLLINS FUNCTION 

After obtaining the double ratios one can try to pa- 
rameterize these measurements in terms of the ColUns 
functions. 



Currently these Gaussian widths have to be taken as 
additional parameters. They are assumed to be univer- 
sal between favored and disfavored fragmentation func- 
tions. In addition, SU(2)/ symmetry for u and d quarks 
is assumed for both types of fragmentation functions: 



A. Assumptions and input 



Assuming a Gaussian dependence on the intrinsic 
transverse momentum 



Di{z,kt) 



cxp(-fc2/(fc?)) 



(27) 



relative to the quark-antiquark axis it is possible to 
solve the convolutions of transverse momenta in the Aq, 
Eqs. dS]),®, asymmetries and to relate them to the A12 
asymmetries, Eqs. Q,®. As the Collins function has to 
obey the following positivity constraint lUF 



Ht{z,h)^<D,{z,kt) 



(28) 



the widths of the Gaussian distributions for the unpo- 
larized fragmentation functions ((fcj)) and the Collins 
fragmentation function ((fej^)) have to differ, since oth- 
erwise the constraint will not hold at sufficiently large 
transverse momenta. 





= Hi 


.U—fTl^ 


= Ht 


,d^7T 




= Hi 




= Hi 


.'U—>7T~ 


Ht^'''{z,kt) 


= Hi 


.U—f-K 


= Hi 






- Hi 


.cl — >n~ 


= Hi 





(29) 
(30) 



and similarly for the unpolarized fragmentation func- 
tions —f Di. The strange quark fragmentation is 
ignored. Under this assumption one can rewrite the dou- 
ble ratio asymmetries entirely in terms of favored and 
disfavored unpolarized and Collins fragmentation func- 
tions. Applying the asymmetry definitions (see Eqn. [5]) 
and integrating over the transverse momenta, the A12 
asymmetries thus become |[26[l (H^ {zi ) is abbreviated as 
Hi and -Di(zi) as Di and equivalently for the Z2 depen- 
dent functions H2 and D2): 
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(32) 



The Aq asymmetries have a similar expression that dif- 
fers by a factor tt/2 tfH]. 



The Gaussian widths of the Collins functions have been 
fixed to be {k^(^)/AP = 2, where M is the pion mass. 
The unpolarized fragmentation functions in the denomi- 
nator are taken either from ["^^KKretzer), ll43h(HKNS). 
ll44](DSS) or [45](KKP) at the scale of = m GeV^ in 
leading order. Since the latter do not contain explicit fa- 
vored and disfavored fragmentation functions they were 
rescaled by (1 + z)/2 for favored fragmentation and 
(1 — z)/2 for disfavored fragmentation functions accord- 
ing to an assumption by 114611 . 



B. Parameterization 

Different parameterizations as a function of the frac- 
tional energy z are possible for the Collins function. The 
simplest case takes Collins functions proportional to z 
times the unpolarized fragmentation functions: 



fr/""(z) = azD{''\z) 



i?f"(z) = bzDf%z) 



(33) 
(34) 



In order to evaluate the sensitivity of such a parameteri- 
zation we calculated the as a function of the two pa- 
rameters a and b. An example is shown in Fig. [24] for 
a combination of Aq^ and Aq'^' asymmetries using the 



KKP unpolarized fragmentation functions 14711. As the 
divided by the number of degrees of freedom (32- 
2) is around 3 a linear description in terms of two pa- 
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rameters seems to be too simple. The parameterization 
does not put stringent constraints on a and b, because 
the minimum of occurs along a diagonal a — b = 
const in the a, b parameter plane. This can be explained 
by the quadratic nature of the measured asymmetries in 
terms of the Collins fragmentation functions. Qualita- 
tively speaking, a very large, favored Collins function 
can be compensated by an almost equally large disfa- 
vored Collins function. While it was hoped that the 
UC data could resolve this ambiguity it turns out that 
the sensitivity is limited and the favored over disfavored 
Collins function ratio H(°'"{z)/Hf''{z) remains inconclu- 
sive. However, if one restricts the ratio of the favored 
Collins function to the favored unpolarized fragmenta- 
tion to be below unity then the opposite signs for the fa- 
vored and disfavored Collins functions suggested by iHl 
can be confirmed. 




FIG. 24: distribution as a function of the parameters a and 
b as specified in the text using the KKP parameterization of the 
unpolarized fragmentation functions for the Ao data. The 
values above Xmax '■= 10 x the minimum were set to 65% of 
Xmax for better visibility. The minimum lies slightly below 100. 
The vertical and horizontal lines correspond to the parameters 
a and 6 of the lowest x^- 



Several models predict the magnitude and the z de- 
pendence of the Collins functions, although favored 
and disfavored fragmentation are in general not dis- 
tinguished. Details can be found in the references 
lllalil,!!!,!^]- While mostly the [1/2] transverse mo- 
mentum moments of the Collins functions are modeled 
these models are consistent with a bare Collins function 
rising with z relative to the unpolarized fragmentation 
function. They are thus also consistent with the rising 
asymmetries presented here. 

In Ref [30] the previously published Belle data was 
used to parameterize the Collins functions. The results 
were then compared to the parameterizations obtained 
from the HERMES |4| and COMPASS H] data assum- 
ing a quark transversity distribution based on the chiral 



quark soliton model. Good agreement between the SIDIS 
and Belle data was found despite the different energy 
scales. In Ref IH] a global fit based on the HERMES, 
COMPASS and the previously published Belle data has 
been performed to extract the first transversity distribu- 
tions. However these parameterizations were still lim- 
ited by the statistics of the previously published data set. 
With the newly obtained results it becomes possible to 
use the Collins functions to further constrain the quark 
transversity distribution. 



VII. SUMMARY AND OUTLOOK 

We have presented a precise measurement of trans- 
verse spin asymmetries, which can be attributed to the 
product of a quark and an antiquark Collins function. 
The statistics has been improved by almost factor of 20 
compared to the previously published results partially 
due to the inclusion of the data taken at the T(4S') res- 
onance. Different combinations of pion pairs exhibit sig- 
nificant, nonzero asymmetries. The systematic uncer- 
tainties in the measurements are understood and eval- 
uated. The results of the measurements show significant 
spin-dependent effects in e+e^ ^ 99 = u, d, s) frag- 
mentation. In addition, assuming a Gaussian transverse 
momentum dependence of the Collins and unpolarized 
fragmentation functions, we studied a possible param- 
eterization of the former Due to the quadratic nature 
of the double ratios the difference between the favored 
and the disfavored Collins functions is still poorly deter- 
mined. However, suggestions based on semi-inclusive 
DIS data, that favored and disfavored Collins functions 
are both large and of opposite sign, are compatible with 
our measurement. A global analysis, which was already 
applied using our previously published data set should 
further constrain both the Collins and the transversity 
functions. 
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TABLE IV: Aq^ and ^fg^ values obtained from fits to pion double ratios as a function of z. The errors shown are statistical and 
systematic. 
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TABLE V: Aq and Ai-f values obtained from fits to pion double ratios as a function of z. The errors shown are statistical and 
systematic. 
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